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Abstract 

A mathematically based model for power MOSFET capacitances of- 
fers both accuracy and speed for circuit simulation applications. 
Expressions for the MOSFET capacitances Cgd(v), C gs (v) and 
Cds( v ) derived from a simplified physical model and 
implemented on the Saber simulator. Data obtained from 60V 
power MOSFETs are used for comparison. 



The objective of this paper is to mathematically describe the capaci- 
tance-voltage characteristics of the interelectrode capacitances in the 
power MOSFET. The internal charge conditions which determine 
the capacitances for different operating states are presented in 
Section 2. In Section 3, experimental three terminal capacitance 
data are compared with simulated results, followed in Section 4 by 
a brief description of how parameters were determined. 



1.0 Introduction 

Numerous MOSFET models have been developed in the past 
decade [8] to replace the generic or standard low voltage MOSFET 
model in SPICE. Most of these models incorporate the generic 
SPICE MOSFET model within a subcircuit of additional passive 
and active elements. However, such subcircuit models are limited 
in flexibility and can be slow to simulate. 

The wide swings of capacitance with terminal voltage typical of 
power MOSFETs make simulations of MOSFET waveforms 
difficult at switching frequencies even as low as 10 kHz. Simple 
power MOSFET models [1-3] assume constant gate-source 
capacitance and a 1 -step or 2-step gate-drain capacitance. A major 
shortcoming of such models is inaccurate modeling of the 
capacitance swings in both the gate-drain and gate-source 
interfaces. A more accurate approach is to utilize a power series 
expression for Cgd [6] but this becomes complex to set up and 
computation- intensive. The most accurate subcircuit model uses a 
special MOS capacitor available in a particular version of SPICE 
[71. A review of existing power MOSFET models [3,8] indicates 
the need for precise expressions of MOSFET capacitances as 
functions of tenninal voltages. 

In this paper the capacitance- voltage eq uatio ns are derived and then 
directly implemented into a power MOSFET model as an alternative 
to the sub-circuit approach. The goal is to develop a fast and simple 
power MOSFET model whose parameters depend on easily 
measured external device properties rather than internal phenomena 
By deriving equations directly from device physics one obtains 
accurate functioning over a wide range of circuit conditions and 
device types. 
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2.0 The Seven Operating States of the Power MOSFET 

All three terminal capacitances Cgd (Vgd)» C gs (Vqs) and 
Cds( v DS) must be modeled together as one system since the 
internal mobile charges which give rise to the capacitances interact 
with each other and are shared among the three terminal 
capacitances. The observed gate capacitance variations are due to 
changes in charge location under the gate in both the n- drain and p- 
body regions. Each of these two regions can be separately in either 
accumulation, depiction or inversion state. When one 
considers all possible combinations of these three states for both 
regions under the gate, one obtains the nine possible operating 
states shown in Table 1. Each state represents a specific 
arrangement of the internal stored charge under the gate of the 
power MOSFET, 

Table 1. The seven operating states of the power MOSFET 
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Closer examination indicates that only seven of the nine states are 
possible for a practical power MOSFET. Six are quasi-static states 
and the seventh, the OFF-4 State, is a dynamic state involving 
charge and discharge of an isolated inversion layer of holes trapped 
at the surface of then- drain regioa The remaining two states: the 
simultaneous accumulation or inversion of mobile carriers in both 
the n- drain and p-body regions, are either impossible or destructive 
to a power MOSFET. The boundaries which separate one quasi- 
static state from another are the following terminal voltages: 

1) V T - Gate-Source Threshold voltage or Turn-on 

voltage of a power MOSFET. 

2) V CBft) - Gate-Body flat-band voltage for the p-body 

region. 

3) VQDfb - Gate-Drain flat-band voltage for the n'drain 

region. 

4 ) v GDI(Vgs) - Gate-Drain voltage at which inversion occurs 

at the n- drain and gate insulator interface. 

The first three are constants, but the gate-drain inversion boundary 
VgdKVgs) shifts with changes in Vcs or V DS bias. At a higher 
drain-source voltage, a more negative (for an n-channel MOSFET) 
gate-source bias voltage is required to induce drain region 
inversion. 

The charge symbols shown in each state of Table 1 are defined as: 
Obi - Body-Inversion charge 
Qbd - Body-Depletion charge 

Qbdc - Body-Depletion constant charge (inversion state) 

Qba - Body-Accumulation charge 

Qdi - Drain-Inversion charge 

Odd - Drain-Depletion charge in OFF-1 , OFF-2 

Qddi - Drain-Depletion charge in OFF-3. OFF-4 (inversion) 

Qda - Drain- Accumulation charge 

Qj - Junction depletion charge of the Body-Drain diode 

These charges are functions of terminal voltages and their equations 
are derived using the set of simplifying assumptions commonly 
known as the Delta-Depletion Approximation [9]: 

a. The depletion region has a uniform charge density containing 
only ionized atoms. 

b. The charges on the gate and interface are contained in 
extremely thin layers. 

c. The charge density is constant along the length of the inversion 
layer. 

Although previous power MOSFET models have assumed that the 
p-body is heavily enough doped to have a constant gate-body 
capacitance, our measurements indicate that the gate p-bodv 
interface contributes a 3:1 variation in C gs with voltage. Thus, the 
complete delta-depletion model must be applied to both the n- drain 
and p-body regions as diagrammed in Figure 1 . 

To individually identify the currents for C g d(v) and C gs (v), the 
total gate charge q G (v) is separated into q G s(v) for current through 
Cgs(v) and qo D (v) for current through C g d(v). Furthermore, the 
highly nonlinear Cds(v) of the rx>dy-drain junction is determined by 
the charge qr>s which depends on both the junction depletion charge 
Qj(Vds), and the inversion charge Qddi in the n- drain region. 
From a circuit perspective, the gate current is the rate of change of 
the total gate charge + 9fipj and me drain currem 

is qGD) For each of the seven practical states in Table 

1. the charge transitions are described next, following the order 



which occurs as a MOSFET tarns off. 
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Figure 1 . Three possible surface charge states for the n-type drain 
region (above) and the p-type body region (below) 

TheON-State of the MOSFET corresponds to accumulation in the 
?5° D ' Wi *lY? D * V GD« and inversion at the p-body 
and gate interface, with V GS £ V T . The internal charge equations 



<Igs = Qbi + Qbdc 
VCD = Qda 
qds-Qj 



(1) 
(2) 
(3) 



Here, Qbdc is a constant charge which does not affect the current 
through C gs (v). It is needed to ensure charge continuity for q GS (v) 
in transition between depletion and inversion. The physical location 
of these charges is shown in Figure 2. 
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Figure 2. MOSFET in the ON-State showing the physical locations 
of the charges Qbi, Qbdc. Qda 

The Active State of the MOSFET occurs when the p-body is 
already in inversion, but a relatively high V D s exists which keeps 
the n- drain in depletion, with V G d < V GD fb. Figure 3 shows that 
during this quasi-static state, charge distribution in the p-body is 
unchanged while q G D becomes: 



<7G£>= Odd 



(4) 
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Figure 3. MOSFET in the Active State showing gate-drain 
depletion Qdd as well as Qbi and Qbdc 

In the OFF-1 State, both the n- drain and p-body regions are 
depleted of mobile carriers, corresponding to the earliest stage of 
the power MOSFETs OFF condition. This state of operation is 
defined by VoBfb < v OS < Vt and Vgd < VcDfb- Here. Figure 4 
shows that the drain depletion region remains unchanged and 



qGS = Qbd 
qcD - Odd 



(5) 
(6) 



the p-body to the source resulting in an increase of Cds(v) due to 
connection of the depletion charge (Qddi) across the drain-source 
electrodes. This quasi-static state is named OFF-3, where Qddi is 
now a function only of Vds. 

qGS = Qba + Qdi + Qddi (8) 

qGD = 0 (9) 

qDS = Qj + Qddi (10) 

Note that inversion in the n" region can occur only after 
accumulation has taken place along the surface of the p-body 
region. 
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Figure 4. MOSFET in the OFF-1 State showing Qbd and Qdd. 

As the gate voltage swings further negative, accumulation of hole 
carriers in the p-body begins to occur when Vgs ^ Voboj. If V D s 
is high, such that Vgd > Vgdi(Vgs)» then the relatively negative 
potential of Qba makes inversion at the n* drain surface impossible 
as shown in Figure 5. We call this the OFF-2 State of operation. 
Here qj>s and qcD are the same as those in OFF-1, but 



qcs = Qba 



(7) 
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Figure 5. MOSFET in the OFF-2 State showing Qba and Qdd- 



If Vds is low and Vgs £ VQBfb. then an inversion layer of holes 
will form at the n- drain surface when Vgd ^ VgdiCVgs)- Once 
this inversion layer forms, the gate-drain is effectively shielded [7] 
and Cgd(v) drops to zero. Figure 6 shows how this shielding 
effect also causes a "short circuit" of the n- inversion region across 



Figure 6. MOSFET in the OFF-3 State showing Qba. Qdi. and 
Qddi- 

In contrast to the other states which are all quasi-static, the OFF-4 
is a dynamic state in which the p-body region is in depletion and 
can trap a layer of inversion charge along the surface of the n- drain 
as shown in Figure 7. This trapped charge will relax to its steady- 
state value in a time constant corresponding to the surface lifetime X. 
The OFF-4 state occurs during MOSFET tum-off and turn-on at 
low drain voltages Vds- With QoKto) the initial inversion charge 
and Qdi the steady-state value, the charge equations are 



qGS - Qbd 



(11) 



qcD = Qddi + QDi(to)ap(- 1 ) + QdAI - «P<- ^ )] (12) 



qos = Qj 



(13) 



Body Depletion Qbd 



Drain Inversion Qdi 




Drain Depletion Qddi 



n- Drain Region 



n+ Drain 



Figure 7. MOSFET in the OFF-4 Stale showing Drain Inversion 
Qddi with Body Depletion Qbd 

No attempt has been made to include the dynamic characteristics of 
the OFF-4 state in the computer model. Instead this state is 
modeled and implemented on the simulator as a quasi-static state, 
assuming the trapped charge instantly relaxes to its steady-state 
value (t = 0 in (12)). The dynamic effect is not considered 
significant enough now to warrant the additional complexity in the 
model. 
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In the Reverse State, lire p-body has yet to achieve inversion and 
remains in depletion. Thus, V GS < V T . However, a low V DS 
allows the n- drain surface to reach accumulation where Vqd £ 
V GDfb- 



QGS = QBD 



(14) 



Here, qcD and qds are me same as Eq.(2) and (3) respectively. 
This state occurs during resonant switching where vds(0 reverses 
polarity while the power MOSFET current ir>s(0 is zero. 

The complete set of power MOSFET charge equations as functions 
of the terminal voltages V GD , V gs . V ds are derived in the 
Appendix. The charge equations are implemented in a three- 
terminal capacitance model, called a template, in the Saber 
simulator. The template reinitializes the simulator whenever the 
simulation crosses a boundary between two of the seven operating 
states of Table 1. The capacitive currents at the three device 
terminals are determined by taking the derivatives of these charce 
expressions. " 



3.0 Experimental Data and Simulation Results 

Capacitance-voltage measurements on a P8NO8 Power MOSFET 
were taken over a complete range of V DS and off-state V GS bias 
conditions. The three interelectrode capacitances Cgd(v), C gs (v) 
and Cds(v) were measured separately using the three terminal 
measurement technique shown in the schematic in Figure 8. 




Figure 8. Circuit for measuring C g d(V GS ,VDs) 

The effects of parallel impedances are eliminated by using the guard 
terminal which functions as the circuit common. The guard 
terminal is electrically different from the ground which is connected 
directly to the instrument chassis. 

With the ( guard properly connected the measurement of the 
unknown capacitance or impedance Z x is unaffected, because Z x 
is calculated using current I x , In Figure 8. the current 12 is 
negligible compared to I x because the internal resistance of the 
ammeter is very low in comparison to the impedance Z2. 

Data on Cds(v) is shown in Figure 9. For a constant V GS bias of - 
5V on the MOSFET gate, a decreasing V DS moves the n- drain 
region from depletion (OFF-2) to inversion (OFF-3). At the 
transition the drain depletion charge Qddi is added to the drain- 
source junction depletion charge Qj resulting in the abrupt increase 
in Cds- The simulated Cd s (v) plot shows the same transition as an 
abrupt jump in the capacitance value at the point where inversion 
occurs. For V GS = 10V and below, both the measured and 
simulated Cds are shown to operate in the OFF-3 state; while at a 
low V GS of OV, the drain-source capacitance operates only in the 
OFF-1 state. The gate-body flatband voltage (V GB n>) has been set 
to -2V for all the simulations The abrupt jump in capacitance value 
as shown for V GS = -5V is a result of the delta-depletion 
approximation assumed for the capacitance model. 




o5> 



Figure 9. Measured Cds (dash) and simulated Cds (solid) versus 
V DS for V GS = OV, -5V, -20V. The arrow shows the transition 
from OFF-2 to OFF-3 as Vds decreases. 




Figure 10. Plots of Cg S (v GS ) for V DS = IV. 10V, 40V. 
Measured (dashed) and simulated (solid). 



Data on Cg S (v) is shown in Figure 10 where for V-r >V G s > 
v GBfb. the n- drain region is depleted of mobile carriers (OFF-1). 
As V G s becomes more negative, inversion takes place along the n~ 
drain interface with the gate, causing the effective area of Cgs and 
therefore, its magnitude, to increase. As V D $ increases, inversion 
in the drain occurs at more negative V G § values. The simulation 
results show that the curvature of the body-depletion capacitance is 
unchanged for all Vps bias in contrast to the varying curvature of 
the measured body-depletion capacitance. For V G s < V G Bfb and 
V GD < V G di(V gs ). the model gives a constant OFF-2 state 
capacitance (flat after the depletion-accumulation boundary at 
V GS = -2V which again is not observed in the measured data. 
These descrepancies are most likely due to the use of the delta- 
depletion model in the simulation. 

Alternatively, if drain inversion occurs prior to Vqs £ V G Bft>» the 
capacitance Cgs(v) moves from OFF-3 to OFF-1 via OFF -4 as Vqs 
is increased. In the model, the OFF-4 state is modeled as a quasi- 
static state which produces an abrupt rise in capacitance during the 
OFF-4 state as shown in Figure 1 1 . Since this is a dynamic state, 
high frequency capacitance measurements will not show the 
capacitance peak in Figure 1 1. However in an actual MOSFET 
switching transient, the charge associated with the large OFF-4 
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capacitance will be transferred as the device moves from the 
dynamic OFF-4 into the quasi-static OFF- 1 or OFF-3 state. For 
this reason, the OFF-4 state is important to retain in the model. 



Table 2. Parameters used in charge equations 




-3 -2.75 -2.3 -3.20 



Figure 1 1. Simulated Cgd(vcs) showing OFF-4 state for V DS =0V 

With V G s = 3V £ V T , the Cgd(vGD) plot goes from the ON-state 
to the Active-state of MOSFET operation as shown in Figure 12. 
The simulated CgdOcD) exhibits a constant capacitance cor- 
responding to Cgdon in the ON-state for low Vds» and then a 
decreasing depletion capacitance in the Active- state for Vds > 
V GS - VoDfb- For Vgs = -3V, an inversion layer at the n" drain 
surface forms during Vgs - V DS ^ VqdiO^Gs) which "shields" the 
gate-drain depletion charge from the gate-drain interface making 
Cgd(v) effectively drop to zero. The simulated Cgd(v) plot again 
shows the transition between OFF-3 and OFF-2 states as an abrupt 
jump in capacitance value as a result of the delta-depletion 
approximation. For Vgs • Vrjs < Vgdi(V*Gs). Cgd(v) depletion 
capacitance decreases as Vds is increased. 
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Figure 12. Plot of Cgd(v D s) for V GS = 3V, -IV, -3V. 
The measured plot is dashed; the simulated solid. 



4.0 Parameter Determination 

The set of MOSFET charge equations in in the Appendix contain 
the nine parameters given in Table 2. This set of parameters are 
required to characterize all three MOSFET interelectrode 
capacitances. 



Five of the seven parameters Cgdon* Cgson, v t, VoBfb and Cjo 
were extracted from the measured capacitance data. Two parameters 
v G Dfb and Vxd w ere determined by curve-fitting the simulated 
capacitance curves with the measured ones. The two built-in 
voltages Vha and V^d were simply assumed to be equal to 0.8 V. 

A more practical and direct way to determine the MOSFET 
parameters is needed. For example, many of the capacitance 
parameters may be determined from the gate-charge plots in 
MOSFET data books. 



5.0 Discussion 

The use of the delta-depletion approximation to model the MOSFET 
capacitances greatly simplifies the modeling of the internal charges 
within the accumulation, depletion, and inversion states. However, 
the complete model with seven operating states becomes complex. 

Nevertheless, the relative accuracy of the capacitance model in 
simulating the real capacitances validates the mathematical equations 
and the assumptions made in creating the model. Furthermore, the 
total CPU time on an H-P 3000 for a complete capacitance sweep 
ranged from 5 to 30 seconds. 

The same delta-depletion concept introduces abrupt changes to Qgs 
and Qgd during inversion of the p-body and the n* drain as demon- 
strated by the abrupt jump in the capacitance value occurring at the 
point of inversion. These charge and capacitance discontinuities are 
a major drawback of this implementation. They can make 
simulation difficult and frequently cause simulation failures. 

Perhaps, this limitation of the delta-depletion model can be 
overcome by a more gradual modeling of the transition between 
depletion and inversion. 

Or alternatively, such a model may be judged too complex for 
practical use in a circuit simulator, and simplification should be 
considered. For example, if the 3:1 variation of Cgs(v) is omitted, 
the complexity of the model can be decreased almost in half. 
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6.0 Conclusions 

A comprehensive charge-based model for the power MOSFET 
capacitances is derived and implemented on the Saber simulator 
Comparison of the simulated capacitances with measured data 
verifies the validity of the capacitance model. 

However, the power MOSFET is a deceptively complex device to 
model accurately, and further work is required to make a practical 
model for circuit simulation purposes. Still needed are: 

1) experimental clarification of the critical performance features 
necessary to include in a power MOSFET model 

2) a simple and direct parameter detenriination method. 

3) an improvement to the modeling of capacitance transitions 
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Appendix: Derivation of the Q(v) Charge Expressions 
The symbols used in this section are: 

Cgdon, Cgson - Constant, ON-State values of the gate-drain 

and gate-source capacitances. 
VY, Vjo - Threshold voltage for the p-body and the n- drain. 
Na , Nf> - doping concentration of both p-type and n-type 

semiconductor. 

XdB, XdD - depth of the depletion in p-body and n~ drain. 
A - area of active region. 

Vna . Vnd - constants (with units of voltage) representing 

fixed physical quantities for the p-type and n- 

type semiconductor. 
VbD - zero-bias built in voltage between the drain and 

surface during inversion. 
VbA - zero-biased built in voltage between the p-body 

and surface during inversion. 
Cjo - junction capacitance of the drain-source interface. 



The charges Qba, Qbi, Qda are simple functions of the terminal 
voltages as can be seen from Figure 1: 



QBI = -CgsoniVGS - Vf) 
Qba = -Cg Son {VGS - V C Bfb) 
Qda = -CgdoniVGD - VGDfb) 



(15) 
(16) 
07) 



All of the above charges go to zero as the voltage approaches the 
boundary value. 

Following the delta-depletion approximation, the n-drain region de- 
pletion charge qDcpieiion of Figure 1 a can be expressed as: 

^Depletion - qNoAXdD = QdD (18) 
This charge is related to the terminal voltages by 



qDeptetion « [qk D (V DS + Vfrp)] 
where kp = 2emNoA 2 . 



(19) 



Next, we want to relate V TD , VcBfb to the physical quantities Nd, 
esi, and Cox of the semiconductor. From the potential equation: 



V T d - V G Dfb = Wax + Vs 



(20) 



where 



Then, Eq. (18), (20). (21) and (22) lead to: 



(22) 



Vtd - V GDJb = (2 V bD y {2V ND j> + V bD (23) 



where 



<* k D 
4C g dor? 



(24) 



Then from Eq. (19) and (24), the depletion charge during inversion 
at the n- drain surface becomes: 



QDDI^Cgdon [*Vnd(VdS + VbD)] 



(25) 



Prior to inversion, the depletion charge is a function of V G d. Thus, 
if we let V DS = V G s - Vqd» 

^Depletion = QdD = -2Cg don V N J^l - - ^y^ 1 ™^ *] 



(26) 



Similarly for the p-body depletion region: 



/ i 

Vt - VGBfb = (2 V bB y (2V NA ? + VhB (27) 



Vna 



(28) 



4Cg SO r? 

QBD = 2C gson V NA [l - (/ + ^GS^Bfb) ^ *] 
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A 



The constant body depletion charge Qbdc during the inversion of 
the p-body is given by: 

Qbdc = 2C gJon v W/t [/ - (# + '] (30) 

During inversion of the n- drain region, the gate charge associated 
with the gate-drain interface is 

<fGD - Cgdoni Vgs ■ V G Djb + Vbo) (3 1 ) 

Since -<JGD = qDcpletion + qinversiom where qjnversion = Qdi, 

Qd! = ~QGD ' QDDI ^ 
= -Cgdor({V GS - V GDfb + V bD ) + [4V ND (VdS + Vft D )] *) 

(32) 

The inversion condition Vgdi( v CS) for the n~ drain region is 
derived by setting Qdi = 0. Then Eq.(29) becomes: 

V G D = V G0 /= V CS + V bD ^^(Vcs - V C Dfb + V bD y (33) 

Eq.(30) shows that at higher drain-source voltages, inversion of the 
n~ drain occurs at increasingly negative gate-source voltages. 

Finally, the charge associated with the drain-source capacitance is 



(7-m) 



CjoVb 
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